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We present density-functional calculations for the interaction energy of monosubstituted 
benzene dimers. Our approach utilizes a recently developed fully nonlocal correlation energy 
functional, which has been applied to the pure benzene dimer and several other systems 
with promising results. The interaction energy as a function of monomer distance was 
calculated for four different substituents in a sandwich and two T-shaped configurations. In 
addition, we considered two methods for dealing with exchange, namely using the revPBE 
generalized gradient functional as well as full Hartree-Fock. Our results are compared with 
other methods, such as M0ller-Plesset and coupled-cluster calculations, thereby establishing 
the usefulness of our approach. Since our density-functional based method is considerably 
faster than other standard methods, it provides a computational inexpensive alternative, 
which is of particular interest for larger systems where standard calculations are too expensive 
or infeasible. 
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I. INTRODUCTION 

With the vast development of computational facilities in the last decade, it is now pos- 
sible to analyze small aromatic systems numerically with high precision. Particular interest 
has emerged for the non-covalent van der Waals interaction, important in many areas and re- 
sponsible for phenomena, such as DNA base pair bonding, protein structure and folding, and 
organic molecule crystallization. While the full DNA structure itself is still out of scope, cal- 
culations for the benzene dimer, which is the simplest prototype of aromatic n-n interaction, 
have been reported. 1-4 A step toward larger systems are calculations for substituent effects 
on the benzene dimer. 5 State-of-the-art tools for such calculations are second-order M0ller- 
Plesset (MP2) perturbation theory and coupled-cluster calculations with singles, doubles, 
and perturbative triples (CCSD(T)). These methods, especially the latter, are regarded as 
highly accurate and reliable. Unfortunately, their application to large systems is time con- 
suming and computationally very demanding. Furthermore, the results might depend on the 
basis set used and elaborate techniques are necessary to estimate the basis set limit. 

Density functional theory (DFT), especially in the form of various generalized gradient 
approximations (GGA) including those with no empirical input, 6,7 have had a degree of 
success in describing dense bulk matter, 8 and isolated molecules. 9 For molecules, this per- 
formance has been further improved by the use of various empirical and hybrid methods. 
However, for van der Waals complexes and sparse matter including many layered struc- 
tures, polymer crystals, organic molecular crystals and the like, such methods either give 
sporadic results or fail completely. In particular, the importance of DFT-based methods 
for calculations of molecular clusters has been minor. The reason for this is that the domi- 
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nant part of the stabilization energy in many cases comes from the dispersion energy, which 
is usually not included correctly in DFT. 10-14 Moreover, functionals of the GGA type of- 
ten do not even predict a binding of molecular dimers. Therefore, considerable effort has 
been devoted to develop DFT based methods that are capable of treating van der Waals 
systems. Among these are methods with empirical damping functions, 15 ' 16 self-consistent 
charge- density functional tight-binding methods, 17 and self-consistent LCAO methods with 
standard pseudopotentials. 18 

Recently, a new approach using a van der Waals density functional (vdW-DF) with a 
nonlocal correlation-energy has been developed. 10 This formalism includes van der Waals 
forces in a seamless fashion and the applicability of this approach was first demonstrated for 
bulk layered systems such as graphite, boron nitride, and molybdenum sulfide. 14 Similarly, 
it has been applied to two-layer dimers of these same materials 19 and the benzene dimer in a 
number of geometries, where it also proved to give a realistic description. 10,20 More recently, 
it has been applied to naphthalene, anthracene, and pyrene dimers 21 ' 22 where for napthalene 
it gives a creditable comparison with recent CCSD(T) results. 23 These larger dimers show 
the power of the method for more sizable systems where CCSD(T) results appear to be 
nonexistent. In addition, there have been promising applications to a molecule physisorbed 
on an infinite surface, 24 as well as to a bulk crystal of polyethylene. 25 

In this paper, we investigate the capability of vdW-DF to describe larger systems. We go 
beyond the pure benzene dimer and use the nonlocal functional to calculate the effect of the 
substituents OH, CH 3 , F, and CN. These each replace an H atom on one of the molecules 
of the benzene dimer, and among other things act as electron donors (OH and CH 3 ) or 
acceptors (F and CN), as discussed by Sinnokrot and Sherrill, 5 who performed MP2 and 
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CCSD(T) calculations. Our results show that vdW-DF, which is much faster than MP2 and 
CCSD(T), gives interaction energies and distances comparable with these standard methods. 
Furthermore, it generates the correct energy ordering of different substituent structures. 
Thus, it is a promising method for larger applications such as DNA base-pair stacking and 
beyond, as well as for various types of extended systems where wave function methods are 
typically impractical. 

This paper is organized as follows: Section II describes the energy functional and the 
details of its evaluation. In Sec. Ill, we give computational details. Our results are analyzed 
in Sec. IV, where we also compare with other numerical methods. We conclude in Sec. V 
and give an outlook for future work. 

II. VAN DER WAALS DENSITY FUNCTIONAL (VDW-DF) 

At the heart of every DFT calculation lies the energy functional E[p], which depends 
on the electron density p. 26 For our calculations we used the following particular functional 
with nonlocal correlation energy, 

E[p] = T s [p] + V pp [p] + J[p] + E x [p] + E^[p] + E^[p\ , (1) 

which was especially designed to include the van der Waals interaction. 10 ' 14 Here, E is 
the total energy functional of the dimer, T s is the single-particle kinetic energy, V pp is the 
ionic pseudopotential functional, and J describes the Coulomb interaction. For the exchange 
functional E x we used the revPBE parameterization of GGA. 27 Furthermore, E\ denotes the 
local part of the correlation, which was evaluated using LDA. 28 The nonlocal part can be 
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written as 

E? L ip] = \ff d 3 rd s r' p(r)cf>(f,r')p(r') , (2) 

where <p(f,f') is a function depending on |r — r'\ and p in the vicinity of r and r'. Details 
about its construction can be found in Ref. [10]. No empirical input was applied. 

For the charge density p, we used the result of a self-consistent DFT calculation utiliz- 
ing the Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional. 7 With this charge 
density we then evaluated the functional from Eq. Differences in interaction energies 
for charge densities resulting from other exchange-correlation functionals turned out to be 
negligible. 

III. COMPUTATIONAL METHODS 

Our self-consistent DFT calculations for obtaining the charge density p were performed 
with the planewave code ABINIT. 29 We used Troullier-Martins pseudopotentials 30 with a 
kinetic energy cutoff of 50 Rydberg. The dimers and monomers were placed in boxes of 
such size that one system and its mirror image is separated by 15 to 20 A of vacuum. 
This distance proved to be sufficient for the charge density to approach zero well before 
the supercell edge, which eliminates any interaction between periodic replica, aside from 
the possibility of electrostatic as discussed below. Even with this conservative box size, 
the calculations of this GGA part go quickly, requiring only around two hours on a single 
processor (e.g. Athlon MP2000+) for each position of the monomers considered. 

Most parts of the energy functional from Eq. (JTJ were also directly evaluated in ABINIT. 
Only the nonlocal correlation part was calculated in a post-process procedure. Thereby, we 
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computed on a real space grid with a spacing equal to the Fourier transform spacing in 
the plane wave calculations, and then we performed a numerical integration as described in 
Ref. [10]. This procedure, as currently coded, requires a computer time comparable to the 
DFT calculation. Although optimization of the code could substantially cut this time, that 
effort is probably unnecessary, as the time for the straightforward GGA calculation of the 
above paragraph will dominate as the system size increases. 

At this point we should mention another promising complementary method (SAPT- 
DFT), 31,32 a density functional based version of symmetry adapted perturbation theory 
(SAPT), which uses time-dependent DFT to calculate the dispersion contribution. This 
complementary procedure appears to have advantages and disadvantages with respect to our 
much simpler one. First, it gives 32 an energy curve closer to CCSD(T) than vdW-DF for the 
benzene dimer 20 in the sandwich geometry and also for rare gas dimers. Furthermore, with 
recent optimization, the computer time required for benzene is (like vdW-DF) comparable 
to that of an ordinary density functional calculation. 33 It appears likely that for systems 
this size or a little larger, this better accuracy may continue. However, the computer time 
of SAPT-DFT scales as the fifth power of the basis size, so that its application will rapidly 
become impractical as the system size increases, as opposed to vdW-DF, which will continue 
to scale as an ordinary DFT calculation. In addition, SAPT-DFT is not a density functional; 
its application necessitates the identification of individual fragments, so it cannot be seamless 
as fragments merge together to become single entities. It is also not obvious how it might be 
applied to a single fragment for which dispersion interaction between different ends or parts 
is important; nor to an extended system that cannot be divided into individual summable 
finite parts. Being a simple density functional, vdW-DF can do all of the above. In summary, 
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it can be applied to calculations of all types and sizes for which an ordinary LDA or GGA 
calculation can. These complementary aspects of SAPT-DFT and vdW-DF render both 
methods important to pursue. 

Because the GGA part of our calculation is performed on a system consisting of replicas 
of the relevant molecular complex separated by 40 Bohr on a simple cubic lattice, one must 
consider the possibility of spurious electrostatic interactions between these replicas showing 
up in our results. Of course, the pure benzene monomer does not have a permanent dipole 
moment. On the other hand, the dipole moments of the substituted benzene range from 
~0.4 debye for toluene to ~4 debye for benzonitrile. The Coulomb energy associated with 
this interaction is easily calculated using a method attributed to Lorentz 34 (see Jackson 35 
for a more easily accessible description). We find that this Coulomb energy ranges from 
~0.05 kcal/mol for benzonitrile and becomes two orders of magnitude smaller for toluene. 
The former number is comparable to the numerical accuracy of our calculations. We cite it 
because it is a hard number that must reasonably be expected to be larger than any induced 
electrostatic interactions, which we will estimate below. This direct interaction obviously 
cancels out of our calculation exactly when the subtraction is made to find the energy of 
interaction between the two parts of each molecular complex. However, the dipole moment 
of the complex in one supercell can induce a changed charge distribution in each of the two 
molecules in an adjacent supercell, thereby changing the electrostatic interaction between 
the latter. This interaction does not cancel out, but can be estimated from the dipole 
moments and polarizabilities involved. We find that it gives a contribution ranging from 
~1CT 4 kcal/mol for the CH 3 substituted dimer to ~1CT 2 kcal/mol for the CN substituted 
one. On the level of the other inaccuracies in our calculation, these errors are completely 
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negligible. 

Of the exchange functional tried earlier, 10 ' 14,19 Zhang and Yang's "revPBE" was the 
closest to Hartree-Fock (HF) for the dimers compared. 10,19 However, it was clear that the 
magnitude of its slope was a little too big in the binding region. This lead to the speculation 
that the equilibrium distances with revPBE exchange, which were also a little too big, would 
be improved if full HF exchange were used instead. This conjecture was verified in Ref. [20], 
where the use of HF exchange instead of revPBE exchange in the functional of Eq. (|T|) 
was fully implemented and applied to the benzene dimer in several different configurations. 
There, the exchange E x was evaluated both with the revPBE exchange functional and HF, 
and the various tradeoffs evaluated and discussed. Here, we do the same for the substituted 
dimers. For HF exchange we used the general ab-initio package GAMESS, 36 using a triple 
zeta valence (TZV) basis set. 

IV. RESULTS AND DISCUSSION 

We calculated interaction energies for three typical configurations of the dimer, 1-4 as de- 
picted in Fig. 1. In the sandwich structure, the two monomers are arranged parallel in such a 
way that the centers of the rings are on top of each other. In the T-shaped structure, the two 
monomers form an angle of 90°, where one monomer is placed right above the center of the 
other ring. We have chosen a configuration in which the substituent points away from the 
other ring [T-shaped(l)] and a similar configuration in which both monomers are rotated 
by 90° [T-shaped(2)]. For our calculations we have considered four different substituted 
molecules, namely phenol, toluene, fluorobenzene, and benzonitrile (whose substituents are 
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OH, CH 3 , F, and CN, respectively). In addition, we also calculated the pure benzene dimer, 
whose "substituent" we refer to as "H". In all cases, we define "distance" as the distance 
from the center of the benzene ring to the center of the aromatic ring containing the sub- 
stituent. All atomic positions for the monomers were optimized during the self-consistent 
DFT calculations, however, for the calculations of interaction energies, the two constituent 
monomers were considered "rigid" and only the intermonomer distance was changed. 

A. Sandwich configuration 

Results for the interaction energy as a function of the distance for the sandwich config- 
uration are plotted in Fig. 2. Each line depicted consists of 20 points for distances from 3 
to 8 A. The data was then fitted with splines to be continuous. As expected, for short dis- 
tances all curves show strong repulsion, whereas for large distances the interaction vanishes. 
Furthermore, all curves show a pronounced minimum around 4A. Values for the interaction 
energy at these minima can be found in Table HI where we also compare with other numerical 
methods. Unfortunately, corresponding experimental values are not available, however, we 
want to mention that some experimental studies of substituent effects in related systems are 
reported, 37 ' 38 but it is difficult to reliably extract correct binding energies. 

It can be seen in Table [Q that the values for the binding distances are approximately 
10% larger than the corresponding MP2 values and about 5% larger than CCSD(T) values. 
Numbers for the interaction energies are in between MP2 and CCSD(T) values, and about 
25% larger compared to CCSD(T). This can be partly explained by the fact that DFT 
calculations underestimate the quadrupole moment by as much as 20%. 39 This makes our 
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interaction energies too large, since not enough repulsion is considered. 

In a simple picture, as pointed out by Sinnokrot and Sherill, 5 one could think of the 
substituents as electron donors (OH, CH 3 ) or electron withdrawing (F, CN). In terms of a 
Hunter-Sanders model, 40 where electrostatics is assumed to determine changes in binding, 
one would expect to find that dimers with F and CN substituents bind more strongly than 
the pure dimer, whereas OH and CH 3 bind more weakly. On the contrary, they found 
increased binding for all the substituted species. A subsequent breakdown 5 using SAPT 
with a modest basis set showed these differences to be distributed among a number of 
different energy components, whose individual values were rather sensitive to the precise 
values of intermonomer distance. We might expect the binding from dispersion to increase 
due to the addition of substituents, as they would be expected to reduce the energy gap for 
electronic particle-hole excitations, and hence increase the polarizability of the substituted 
monomer. However, SAPT calculations showed this effect to be too small to alone explain 
the differences. Nevertheless, it was clear that dispersion was by far the largest individual 
attractive component of the energy, and the necessity to be able to calculate this with little 
variation in fractional error for each substituted species is crucial to obtain a correct overall 
picture. 

In the light of the above, it is an important challenge to the vdW-DF density functional to 
see if it can reproduce the changes that occur upon substitution, and hence it is interesting 
to look at our predicted interaction energies relative to those of the pure benzene dimer. 
Numbers for these values at the corresponding minimum are collected in Table |H] If we also 
normalize the distance such that it is measured relative to the benzene binding distance, the 
plot in the right panel of Fig. 2 is obtained. Here, we also added the CCSD(T) minimum 
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points so that a direct comparison is possible. It can be seen clearly that the DFT method is 
in very good agreement with the CCSD(T) calculations. Similarly, we predict a bond length 
contraction of ~0.1 A for all the substituents, in agreement with the CCSD(T) calculations. 
In general, while all interaction energies depend somewhat upon the specific details of the 
calculation, the relative energies for the different substituents are almost independent thereof. 

B. T-shaped(l) and T-shaped(2) configuration 

Results for interaction energies of T-shaped(l) and T-shaped(2) dimers are plotted in 
Figs. 3 and 4. All curves show pronounced minima around 5.25 A. Values for E- mt at the 
minimum points and relative to the benzene dimer can again be found in Tables H] and El 
respectively. It can be seen immediately that the effect of substitution in the T-shaped 
configurations is much weaker. The difference in E int for the strongest and weakest binding 
substituents is only approximately half of what it is in sandwich configuration. 

As mentioned before, the quadrupole-quadrupole interaction within DFT is underesti- 
mated. For sandwich structures, this leads to a binding that is too strong, since not enough 
repulsion is considered. For T-shaped configurations, however, the quadrupole-quadrupole 
interaction constitutes an attraction, since the monomers form an angle of 90°. Here, the 
underestimation will lead to a binding which is too weak. The effect will be smaller though, 
because the T-shaped binding distances are larger than the sandwich ones. As a result, the 
vdW-DF interaction energies now show a weaker binding compared to MP2 and CCSD(T). 
The binding order, on the other hand, is in perfect agreement. In particular, the isoenergetic 
levels of H and OH in T-shaped(l) and T-shaped(2) configuration are well reproduced. 
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As for the sandwich structure, we calculated the interaction energies in reference to the 
benzene dimer and its binding distance. The results are depicted in the right graphs of Figs. 3 
and 4. Again, it can be seen that the DFT method is generally in quite good agreement 
with the CCSD(T) calculations. 

C. Hartree-Fock exchange 

As pointed out in Ref. [20], using HF exchange instead of revPBE while evaluating 
the functional of Eq. influences binding distances and interaction energies. While the 
exchange energy might alter the shape of the potential curve to a smaller or greater extent, 
it is the nonlocal correlation energy functional that allows to treat van der Waals systems 
within a DFT approach. As we will see below, the binding order of different substituents is 
actually independent of the exchange functional used and thus, it is governed by the nonlocal 
correlation energy and other components of the total energy. From this point of view, it seems 
conceivable to utilize even other exchange functionals than revPBE or HF. However, there 
is one important point to keep in mind: Within a given covalently bonded monomer the 
full HF interaction can be long ranged, because the Coulomb interaction can extend from 
one end of a molecular orbital to another, whereas the proper inclusion of correlation will 
shield this interaction, yielding an effectively short-ranged Coulomb interaction (aside from 
van der Waals effects). In ordinary DFT (either LDA or GGA) this is handled by making 
both exchange and correlation short-ranged, thus mixing their definitions in this respect. 

We repeated all calculations from above using HF exchange instead of revPBE. The 
results for all configurations are plotted in Fig. 5. In general, the effect of using HF exchange 
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instead of revPBE is a stronger binding at shorter distances. Values for the interaction energy 
at the minimum points can be found in Tables ITlandlTTl where we refer to the DFT calculations 
with HF exchange as vdW-DF(HF). Remarkably, the relative order of interaction energies 
for different substituents is the same with HF and revPBE exchange. It can also be seen that 
the effect of using HF exchange instead of revPBE is stronger in the T-shaped configuration 
than in the sandwich structure. However, it is clear that the energy differences between the 
benzene dimer and the four substituted ones are not in such good agreement with CCSD(T) 
as the results using revPBE were. As the substituted monomers are significantly longer and 
presumably more easily polarized than benzene alone, one can speculate that this poorer 
agreement is due to the lack of long range shielding of the Coulomb potential in the HF 
approximation, as discussed in the above paragraph. Alternatively, it could be due to the 
lack of self-consistency of simply substituting the self-consistent HF exchange energy for 
the self-consistent GGA one, which is a related problem. In any case, the larger issue of 
the consistency of various exchange approximations with the correlation method used is an 
important one, and a subject of further study. 

However, with respect to the grosser effects of HF exchange, we find results similar to 
Ref. [20], in particular an improvement in the binding distances. The vdW-DF(HF) values 
are now within ~1% of the CCSD(T) values. On the other hand, the interaction energies 
now show a stronger binding than the vdW-DF results. The use of HF exchange also corrects 
the relation of interaction energies between T-shaped and sandwich structures: While vdW- 
DF always predicts the sandwich structure to be lower in energy, vdW-DF(HF) reverses this 
order and gives stronger binding for all T-shaped structures, in agreement with CCSD(T) 
calculations. 5 
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V. CONCLUSIONS 

We have presented DFT-based calculations for systems that exhibit van der Waals type 
bonding. In particular, we investigated substituent effects on benzene dimers for several 
different substituents. Our method gives useful approximations for interaction energies and 
binding distances, when compared with state-of-the-art wave function methods, which are 
computationally much more demanding. In addition, the changes produced by the sub- 
stituents as given by our vdW-DF method were quite close indeed to those predicted by 
CCSD(T). Our DFT based method offers a computationally inexpensive alternative that 
can speed up calculations drastically, and in addition make applications to extended systems 
possible. To be specific, our calculations require only slightly more computer time for large 
systems than the GGA type of DFT calculations that they replace. With this in mind, it 
is a promising candidate for larger applications, such as DNA base-pair stacking, which is 
already work in progress. The method is also applicable to extended systems where wave 
function methods are infeasible. 
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TABLE I: Comparison of interaction energies E int and binding distances R for all sub- 
stituents, calculated with different methods. MP2 and CCSD(T) values are taken from 
Ref. [5], where calculations were performed with an aug-cc-pVTZ basis set. 





sandwich 


T-shaped(l) 


T-shaped(2) 


method 


R [A] E [kcal/mol] 


R [A] E [kcal/mol] 


R [A] E [kcal/mol] 



H 


MP2 


3.70 


-3.26 


4.89 


-3.46 


4.89 


-3.46 




estd. CCSD(T) 


3.90 


-1.80 


4.99 


-2.62 


4.99 


-2.62 




vdW-DF 


4.10 


-2.37 


5.28 


-2.28 


5.28 


-2.28 




vdW-DF(HF) 


3.83 


-2.93 


4.90 


-4.40 


4.90 


-4.40 


OH 


MP2 


3.60 


-3.75 


4.90 


-3.42 


4.90 


-3.52 




estd. CCSD(T) 


3.80 


-2.17 


5.00 


-2.58 


5.00 


-2.67 




vdW-DF 


4.03 


-2.79 


5.27 


-2.28 


5.26 


-2.41 




vdW-DF(HF) 


3.74 


-3.85 


4.88 


-4.45 


4.90 


-4.39 


CH 3 


MP2 


3.65 


-3.96 


4.90 


-3.39 


4.80 


-3.89 




estd. CCSD(T) 


3.80 


-2.27 


5.00 


-2.55 


5.00 


-2.95 




vdW-DF 


4.04 


-2.95 


5.25 


-2.27 


5.24 


-2.62 




vdW-DF(HF) 


3.75 


-3.81 


4.87 


-4.26 


4.88 


-4.82 


F 


MP2 


3.70 


-3.81 


4.90 


-3.61 


4.90 


-3.17 




estd. CCSD(T) 


3.80 


-2.29 


5.00 


-2.77 


5.00 


-2.38 




vdW-DF 


4.03 


-2.85 


5.28 


-2.46 


5.29 


-2.23 




vdW-DF(HF) 


3.75 


-3.98 


4.89 


-4.81 


4.91 


-4.05 


CN 


MP2 


3.60 


-4.86 


4.80 


-4.11 


4.90 


-3.08 




estd. CCSD(T) 


3.80 


-3.05 


4.90 


-3.25 


5.00 


-2.20 




vdW-DF 


3.97 


-3.68 


5.22 


-2.86 


5.31 


-2.01 




vdW-DF(HF) 


3.70 


-5.16 


4.86 


-5.37 


4.92 


-3.85 
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TABLE II: Interaction energies relative to the benzene dimer 
pVTZ and CCSD(T)/aug-cc-pVTZ values are taken from Ref. 


in [kcal/mol]. MP2/aug-cc- 
[5]. 




OH 


CH 3 


F 


CN 






sandwich structure 






MP2 


-0.49 


-0.70 


-0.55 


-1.60 


estd. CCSD(T) 


-0.37 


-0.47 


-0.49 


-1.25 


vdW-DF 


-0.42 


-0.58 


-0.48 


-1.31 


vdW-DF(HF) 


-0.92 


-0.88 


-1.05 


-2.23 




r 


r-shaped(l) structure 






MP2 


0.04 


0.07 


-0.15 


-0.65 


estd. CCSD(T) 


0.04 


0.07 


-0.15 


-0.63 


vdW-DF 


0.00 


0.01 


-0.18 


-0.58 


vdW-DF(HF) 


-0.05 


0.14 


-0.41 


-0.97 




r 


r-shaped(2) structure 






MP2 


-0.06 


-0.44 


0.29 


0.39 


estd. CCSD(T) 


-0.05 


-0.33 


0.24 


0.42 


vdW-DF 


-0.13 


-0.34 


0.05 


0.27 


vdW-DF(HF) 


0.01 


-0.42 


0.35 


0.55 
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FIG. 1. Models of the dimer, where we have used a CH 3 substituent as example 
in a) sandwich configuration, b) T-shaped(l) configuration, and c) T-shaped(2) 
configuration. The distance is measured between the centers of the rings. 

FIG. 2. Results for the interaction energy E- mt as a function of the monomer 
separation in sandwich configuration. In the left panel, the absolute energy 
values are plotted over the absolute distance. In the right panel, the energy 
and distance are plotted with respect to the minimum energy and distance of 
the benzene dimer. In addition, the CCSD(T) minimum points are depicted as 
corresponding larger symbols. 

FIG. 3. Results for the interaction energy E- mt as a function of the monomer 
separation in T-shaped(l) configuration, (left) The absolute energy and distance 
are plotted, (right) Energy and distance are given in reference to the minimum 
energy and distance of the benzene dimer. In addition, the CCSD(T) results are 
depicted. 
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FIG. 4. Results for the interaction energy E int as a function of monomer sep- 
aration in T-shaped(2) configuration, (left) The absolute energy and distance 
are plotted, (right) Energy and distance are given in reference to the minimum 
energy and distance of the benzene dimer. The CCSD(T) minimum points are 
also depicted as corresponding symbols. 



FIG. 5. Interaction energy E- mt as a function of the distance for sandwich, T- 
shaped(l), and T-shaped(2) configuration, where Hartree-Fock exchange was 
used. 
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T. Thonhauser et al, Fig. 1 




T. Thonhauser et al, Fig. 2 
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